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ABSTRACT 

We use high-resolution adaptive optics (AO) imaging on the Keck II telescope to study the grav- 
itational lens B0128+437 in unprecedented detail, allowing us to resolve individual lensed quasar 
components and, for the first time, detect and measure properties of the lensing galaxy. B0128+437 is 
a small separation lens with known flux-ratio and astrometric anomalies. We discuss possible causes 
for these anomalies, including the presence of substructure in the lensing galaxy, propagation effects 
due to dust and a turbulent interstellar medium, and gravitational microlensing. This work demon- 
strates that AO will be an essential tool for studying the many new small-separation lenses expected 
from future surveys. 

Subject headings: galaxies: high-redshift — galaxies: individual (CLASS B0128+437) — gravitational 
lensing: strong — infrared: galaxies 



1. INTRODUCTION 

Galaxies acting as strong gravitational lenses can pro- 
vide important insights into the distribution of matter 
on small scales, as the brightness and position of the 
multiple lensed images are directly related to the gravi- 
tational potential of the foreground galaxy. Many lenses 
are well modeled by a smooth potential. For other lenses, 
however, a simple mass model is incapable of accurately 
predicting the positions of the multiple images of the 
background object (astrometric anomalies), or the rela- 
tive brightness between image pairs (flux-ratio anoma- 
lies). The existence of these astrometric and flux-ratio 
anomalies in lenses has been used to study the pres- 
ence of galacti c substructure (jDalal fc Kochanekl 120021 : 
Ku et aLll2009| ) - predicted by numerica l simulations of 
the ACDM theory of struc ture formation (jDiemand et al.l 
l2008tlSp"ringel et aT]|2008[ ) - as well as propagation effects 
such as dust extinction and scatter broadening. 

In this letter we discuss one such an omalous lens, 
CLASS B0128+437 (IPhillips et al.l [200(1 . First ob- 
served as part of t he Cosmic Lens All -Sky Survey (e.g., 
iMvers et al.l[200l iBrowne et aDl200l . B0128+437 is a 
compact four-image system with a maximum image sep- 
aration of 0.54". Very Lo ng Baseline Interferom etry 
(VLBI) imaging presented bv lBiggs etHI (|2004; B04fl re- 
vealed three distinct subcomponents in each of the lensed 
source images. Lens modeling in B04 could not suffi- 
ciently reproduce either the positions of these subcom- 
ponents in images B and C or the total flux of image B, 
suggesting that some form of substructure was present 
i n the lens system. 

iMcKean et al.l ((2004) determined a redshift of z s = 
3.124 for the source galaxy, and an emission line not as- 
sociated with the source was posited to be either Ha 
(implying a lens redshift z; = 0.218), H/3 (zi — 0.645), or 
[O II] (zi = 1.145). Their analysis suggested that [O II] 
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was the most likely candidate; if confirmed, this would 
indicate that B0128+437 is the most distant known grav- 
itational lens. 

Given the small image separation and a possible [O II] 
emission line, the lensing galaxy in the B0128+437 sys- 
tem may be a late-type galaxy. Ground based images 
in the infrared (IR) from UKIRT (B04) have not been 
able to resolve th e system, an d space based images in 
the IR (NICMOS; [Biggi|200l Figure Q] left) and opti- 
cal (WFPC2) have had neither the resolution nor, in the 
case of WFPC2, sufficient sensitivity to investigate the 
properties of the system. 

In this Letter we present a different method of ob- 
serving B0128+437: ground-based near-IR imaging cou- 
pled with laser guide-star adaptive optics (LGS AO). 
While AO i maging has been used on gravitational lenses 
in the past dCrampton et al.lll998t iMarshall et al.1120071 
IMcKean et al.l 120071 : lAuger et all 120081 ). the systems 
studied have been high-mass lenses producing image sep- 
arations of -2". With B0128+437, we present the first 
results of using AO with a small-separation gravitational 
lens, demonstrating that AO imaging is capable of prob- 
ing the details of low-mass systems. 

Throughout this Letter, we assume a concordance cos- 
mological model, with Q,m = 0.3, ft\ = 0.7, and a Hub- 
ble parameter Ho = 70 ft.70 km s _1 Mpc -1 . All magni- 
tudes presented are AB magnitudes. 

2. OBSERVATIONS 

We observed the B0128+437 system using the Near In- 
frared Camera 2 (NIRC2; K. Matthews et al., in prepa- 
ration) on the Keck II telescope along with the the 
LGS AO system, under photometric conditions, on UT 
2009 September 12. We used the narrow camera (10" x 
10" field of view, 0.01"/pixel scale) to properly sample 
the diffraction-limited core of the point-spread function 
(PSF). The data were obtained using the Kp filter and 
were dithered to mitigate the effects of bad pixels, cos- 
mic ray strikes, and detector persistence. Each exposure 
was 180 seconds to stay within the linear regime of the 
detector and to minimize the variation in sky brightness 
for flat-fielding. After visually inspecting cleaned indi- 
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Table 1 

Lens System Parameters 



Component 



Kp 



Aa(") a 



A<5(") a 



Image A 
Image B 
Image C 
Image D 

Lens galaxy (total) 
Lens galaxy (i?Ein) b 



21.55 ± 0.03 
23.49 ± 0.22 
22.49 dz 0.04 
22.87 ± 0.12 
20.99 ± 0.08 
22.15 ± 0.10 



= 0.0 
0.099 
0.521 
0.109 
0.217 
0.217 



= 0.0 
0.095 
-0.170 
-0.260 
-0.104 
-0.104 



a Uncertainties are ~0.001 for B/C/D and ~0.01 for the galaxy. 
b Aperture magnitude with R = i?Ein = 0.24". 

vidual exposures and rejecting frames where the image 
appeared strongly distorted due to loss of LGS lock or 
problems with the AO wave front sensors, we were left 
with 30 frames and a total exposure time of 5400 s. 

The data were reduced using the Center for Adaptive 
Optics Treasure Survey (CATS NIRC2 pipeline, mod- 
ified to properly correct for the geometric distortion of 
the narrow camergQ. The reduced i mage frames were 
then stacked using the Drizzle routine (|Fruchter fc Hookl 
2002) in IRAF, using sub-pixel offsets determined via 
cross-correlation. 

The results of the reduction can be seen in the middle 
panel of Figure [TJ We see the four lensed quasar im- 
ages clearly, along with the lensing galaxy and a partial 
Einstein ring from the quasar host galaxy. This is a sig- 
nificant improvement over the 5440 s NICMOS F160W 
image. 

3. MODELING MASS AND LIGHT 

We are primarily interested in the light of the lensing 
galaxy and any substructure that might be associated 
with it and we therefore treat the lensed quasar and its 
host galaxy (the Einstein ring) as nuisances. We model 
the flux in the ima ge as various components: a Sersic 
profile (|Sersid ll963) for the lensing galaxy, PSF models 
for the quasi-stellar object (QSO) images, and a lensed 
Sersic profile for the QSO host galaxy. We then marginal- 
ize over the parameters of the source components to in- 
fer the properties of the lensing galaxy. Note that we 
do not treat the QSO images as being lensed; this is 
because we suspect that simple lens models will not be 
sufficient to reproduce the observed fluxes, complicating 
our attempts to model the light. We do not use the QSO 
image positions to constrain the lens model, because we 
are able to obtain a robust subtraction of the Einstein 
ring by using the data in the ring itself and assuming the 
Singular Isothermal Ellipsoid (SIE) + external shear pa- 
rameters from B04 as a starting point for the lens model; 
the final model is consistent with the B04 model, given 
the uncertainties. It is interesting that we are able to 
sufficiently remove the Einstein ring flux with a simple 
model, as it implies that the known anomalies present in 
the B0128+437 system are not the result of incorrect as- 
sumptions about the global mass model, but rather due 
to some localized discrepancy (§4.2). 

The PSF model is determined from the data by fit- 
ting a model of three Gaussian components to the QSO 
images; one component represents the diffraction-limited 

3 http:/ /irlab. astro. ucla.edu/cats/index.shtml 

4 http:/ /www. astro. caltech.edu/~pbc/AO/distortion.pdf 



core, another represents the seeing-limited diffuse PSF, 
and the third encodes structure in the PSF, as seen in 
the C and D components in Figure [1] (middle). The full 
model (Sersic lens, four QSO PSFs, and lensed QSO 
host) is fit to the data iteratively, using a Levenberg- 
Marquardt optimization algorithm to find the minimum 
in the parameter space. A set of Markov Chain Monte 
Carlo simulations are run to probe the parameter space 
around this minimum and obtain uncertainty estimates 
for the quasar image fluxes, lens flux, and lens position. 
The results of the modeling are shown in Figure Q] (right) 
and the best-fit parameters for the lensing galaxy and 
QSO positions and fluxes are tabulated in Table [TJ 

4. DISCUSSION 

The improved image resolution provided by the AO 
data allows us to (1) for the first time observationally 
detect the lensing galaxy and estimate its light profile, 
(2) measure flux ratios between pairs of lensed images 
without worrying about blending effects, and (3) search 
for luminous substructure. 

4.1. The Lens Galaxy 

In the rightmost panel of Figure [JJ we clearly see the 
BO 128+437 lens galaxy. We measure a total magnitude 
of Kp = 20.99 ± 0.08. Using the B04 model Einstein 
radius i?Em = 0.24", we further measure an Einstein 
radius aperture magnitude Kp.Ein = 22.15±0.10. We are 
unable to obtain robust constraints on either the effective 
radius or the Sersic index, although we note that in all 
iterations of our modeling, the Sersic index lies between 
n = 1 and n = 2.5, which suggests either a lenticular or 
late-type g alaxy morphology. 

We use Maraston (2005; M05) instantaneous burst 
stellar population models with solar metallicity to in- 
fer the age of the lens galaxy stellar population. We use 
models with predominantly stellar lensing mass (i.e., a 
stellar fraction /* ~ 1) within the Einstein radius, which 
are cons istent with th e resu lts for low-mass SO galax- 
ies from iTortora et al.1 (120091) and with previous results 
from lenses (jAuger et al.ll2009l ) after extrapolating to the 
mass of B0128+437. We then infer the observed Kp-band 
magnitude for each of the candidate lens redshifts from 
the M05 models for a range of stellar population ages 
and compare this in a x 2 sense with the observed magni- 
tude within the Einstein radius (Figure [2]). For redshift 
z\ = 0.218, no acceptable solution - i.e. one in which 
the stellar age was less than the age of the universe - 
is found unless we assume an unusually low (/* = 0.3) 
stellar fraction. This analysis does not strongly favor ei- 
ther of the other two proposed lens redshifts, z — 1.145 
and z = 0.645, both of which give stellar ages that are 
less than the age of the universe for reasonable values 
of /*. However, after combining this with the previ- 
ous photometric and spectroscopic data, we continue to 
slightly favor the higher redshift z = 1.145. The inferred 
age for z = 1.145 is —3.2 Gyr (Table E} and the rest- 
frame K-band luminosity inferred from the M05 model 
is log 10 (LK Ein) = 10.61 L Q (the age and luminosity for 
/* = 0.9 are 2.7 Gyr and log 10 _(L K ,Ein) = 10.67 L Q ). This 
results in a mass-to-light ratio M/LK,Ein = 0.68 which 
is, by construction, consistent with a stellar population 
that formed ~3 Gyr ago. 
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Figure 1. Left: NICMOS F160W image of B0128+437. The resolution of HST in this band makes it difficult to resolve individual 
components. Middle: AO image of B0128+437. The four lensed quasar images [A-D] are clearly seen and easily resolved, along with a 
partial Einstein ring from the quasar host galaxy. The 3-cr 2.3 GHz contours of B04 are overlayed in white. Additionally, the lensing galaxy 
can be seen lying close to A, B, and D in the image plane. Right: Flux remaining in the B0128+437 system after subtracting the source 
quasar and host galaxy light. The lens galaxy is now clearly visible. Note that the image stretch in this panel is different from that of the 
middle panel. 



Table 2 

Redshift Dependent Lens Galaxy Properties 





Stellar Age a 


'°SlO ^K.lens 


l°gl0 ^K,Ein 


log 10 M Ein b 


M/L K , Ein 






(hro Gyr) 


(h 70 Lq) 


(h 70 Lq) 


(hro Mq) 


{h 70 Mq/Lq) 


(km s _1 ) 


0.218 


10.8 


9.53 ± 0.03 


9.07 ± 0.04 


9.78 


5.13 ± 0.04 


98 


0.645 


8.5 


10.55 ± 0.04 


10.09 ± 0.04 


10.20 


1.29 ± 0.04 


114 


1.145 


3.2 


11.07 ± 0.05 


10.61 ± 0.05 


10.44 


0.68 ± 0.05 


138 



a Measured assuming a stellar fraction /„ = 1.0 for z\ = {0.645, 1.145} 
b Determined using the B04 mass model with -F?Ein = 0.24". 




Age (Gyr) 

Figure 2. Inference on the stellar age of the lensing galaxy by 
comparing the observed magnitude (Kp = 22.15) with M05 mod- 
els of varying ages. The solid (dotted) curves represent a stellar 
fraction /* = 1 (/„ = 0.9), while the vertical dashed lines denote 
the age of the universe at the assumed lens redshift. The two sets of 
curves represent lens redshifts of z; = 1.145 (black) and 2j = 0.645 
(red); the lens redshift z; = 0.218 leads to stellar population ages 
older than the age of the universe at that redshift for /„ > 0.3. 

4.2. Flux Ratios 

Our models give accurate fluxes for each image sepa- 
rately, allowing us to calculate flux ratios. We present 
these flux ratios, along with the radio flux ratios of B04, 
in Table [3] While there is significant variation between 
the flux ratios in different bands, one feature that per- 



and /* = 0.3 for z t = 0.218. 



Table 3 

Multi-band Flux-ratio Data 



Flux Ratio 


Kp 


X 


C 


S 


L 


f B /f.A 


0.17 ± 0.03 




0.56 


0.49 




Scl Sa 


0.42 ± 0.01 




0.49 


0.34 




fD/fA 


0.30 ± 0.03 




0.47 


0.47 




fc/fD 


1.43 ± 0.16 


0.67 


1.04 


0.72 


0.75 



Flux ratios between image pairs. We present Kp flux ratios, along 
with previously measured flux ratios in the X (3.5 cm), C (6 cm), 
S (13 cm), and L (21 cm) radio bands. In the X and L bands, only 
images C and D are sufficiently deblended to measure accurate flux 
ratios. Due to the extreme blending of flux between images, flux 
ratios measured in the F160W NICMOS image are unreliable, and 
thus not included. All radio data are taken from B04, with errors 
< 5%. 

sists in both the radio and near-IR regimes is the un- 
usual value of /b/Za- The close proximity of A and B 
suggests that they should h ave nearly the same flux and 
a negligible time delay (e.g. JKeeton et aLll2005D : the dis- 
crepancy between what is seen in the data and what is 
expected from smooth mass models therefore cannot be 
due to intrinsic variability in the source. We therefore 
discuss several other phenomena that may contribute to 
the flux-ratio anomaly. 

4.2.1. Substructure 

Substructure, in the form of satellite galaxies and 
dark matter clumps associated with the lensing galaxy, 
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could explain the anomalous /b/Ja, especially since 
some lenses with flux-ratio anomalies are known to 
have luminous satellites whose presence ca n account 
for t he anomalies (e. g., MG041 4 +053 4: iRos et al.1 
20001: MG201 6+112: Mo re et all [200l B2045+265: 
McKean et al.l 12007}) . This scenario is appealing because 
VLBI imaging of BOf 28+437 shows astrometric anoma- 
lies of the quasar image positions. The presence of sub- 
structure perturbs the gravitational potential and con- 
sequently alters the position and the brightness of the 
lensed images; substructure could therefore be responsi- 
ble for both anomalies. 

However, since lensing is an achromatic effect, the mea- 
sured flux ratios in the radio regime should be the same 
as those measured in the near-IR. Instead, the near-IR 
Ib I j A is significantly smaller than its radio wavelength 
counterparts (Table [3]). Therefore, it seems unlikely that 
substructure alone is responsible for the anomaly. 

Furthermore, after subtracting both the source and 
lens flux, we find no obvious luminous substructure in 
the vicinity of B0I28+437; the flux that remains in the 
right panel of Figure Q] is most likely an artifact of model 
subtraction. However, the fact that we do not directly 
detect luminous substructure may not be surprising. B04 
estimate that a perturbing mass as low as 10 6 M could 
account for the milliarcsecond-scale astrometric discrep- 
ancies observed in B0128+437. If we assume reasonable 
M/L ratios for such an object, it falls well below our 
detection threshold. 

We determine a limiting magnitude for our imaging 
by populating the image with randomly placed simu- 
lated point sources covering a range of magnitudes from 
Kp = 23.5 to Kp = 28.5. We then analyz e these images 
with SExtractor (jBertin fc Arnoutsl 119961 ) and find that 
we can repeatedly recover sources as faint as Kp = 25.5, 
which we take to be our detection limit. This corresponds 
to a rest-frame V band luminosity My = —17, compa- 
rable to the Magellanic Clouds. We note, however, that 
the detection limit is significantly brighter very close to 
the lensed images due to the increased shot noise from 
the quasars. 

4.2.2. Propagation Effects 

Propagation effects can also account for the Jb/Ja 
anomaly, with scatter broadening by the lensing galaxy 
interstellar medium affecting the radio fluxes and dust 
extinction affecting the near-IR fluxes, giving rise to 
the chromatic nature of the flux-ratio anomaly. If dust 
were solely responsible for the near-IR anomaly, image B 
would need to be dimmed by ^2.0 magnitudes of extinc- 
tion (assuming as a low er limit that image A is not be- 
ing affected). Using the lCalzetti et all ()2000f ) reddening 
law for rest- frame J-band, this corresponds to an E(B-V) 
value of 1.01, equivalent to ^4.1 magnitudes of rest-frame 
V-band extinction. While this is large, it is not unphys- 
ical: studies of the local group have revealed pockets of 
dust in nearby giant molecular clouds (GMCs) that cause 
up to 10 V-band magnitudes of extinction (jCambresyl 
[19991) . 

In the radio, scatter broadening by a GMC may ex- 
plain why image B shows a decrease in surface bright- 
ness and why its subcomponents seen in 5 GHz VLBI 
imaging are much broader than those of other im- 
ages. This scenario has been postulated for the grav- 



itational lens B0218+357, where a GMC is thought to 
be responsible for the broadening /dimming of im age A 
(|Grundahl fc Hiorthl[l995l IZeieer fe Darling|l2010D . 

However, while this explains the flux-ratio anomaly, 
known GMCs are not massive enough to significantly 
alter the lensing potential and cause the astrometric 
anomalies shown by images B and C. Therefore, it is 
unlikely that propagation effects alone can explain all of 
the observed properties. 

4.2.3. Microlensmg 

A third phenomenon that could explain the fs / /a 
flux-ratio anomaly is extragalactic microlensing, where 
small gravitational distortions due to star s in the lens- 
ing g alaxy affect the fluxes of images (e.g.. iWambsganssl 
2006). While the size of the radio-emitting regions of 
active galaxies are often too large to be significantly af- 
fected by the lensing cross section of individual stars, mi- 
crolensing of the accretion disk (the source of the near-IR 
emission in B0128+437) can substantially change the ob- 
served flux of a given image. Microlensing could, there- 
fore, explain the discrepant flux ratios between the radio 
and near-IR data. 

However, while it is possible to observe some mi- 
crolensing in radio bands, the flux discrepancies mea- 
sured for other lens systems are much less dramatic 
than the radio /b/Ia val ues of B0128+4 37. Stud- 
ies involving B1608+656: dFassnacht et al l l200l and 
B1600+434: (|Koopmans fe de Bruvnl l2000f ) show only 
small (~3%) microlensing radio flux-ratio anomalies. 
Furthermore, based on the movement of stars relative 
to the lensed images, the microlensing signal should be 
ti me- variable. However, rad io monitoring of B0128+437 
by iKoopmans et al.l (|2003[ ) has shown a nearly time- 
independent flux ratio. Finally, microlensing cannot ex- 
plain the observed astrometric anomalies. Thus, it would 
be difficult for microlensing alone to fully account for the 
observed properties of B0128+437. 

4.2.4. Multiple Processes 

After comparing our data to the existing radio data, 
we find that no single phenomenon seems to completely 
explain the unusual appearance of B0128+437. How- 
ever, since each process is able to explain some aspect 
of the anomaly, the full effect is likely due to a combi- 
nation of these mechanisms. Considering the probable 
late-type morphology of the lensing galaxy, this multi- 
process scenario could, for example, be the result of a spi- 
ral arm. The overdensity of an arm positioned near image 
B could be sufficient to cause the astrometric anomaly 
(|Mao fc Schneidenfl998l) . and a GMC that is associated 
with the enhanced star formation in the spiral arm could 
produce the wavelength-dependent flux-ratio anomaly. 

However, additional data and modeling are needed be- 
fore a definitive explanation for the observed properties 
of this system can be made. Obtaining high-resolution, 
multi-band imaging in the optical regime will allow for 
better constraints on extinction. Rest-frame U- and B- 
band data may also provide visual confirmation of spi- 
ral arms, though given the faintness of the source, deep 
imaging using either AO or the JWST will be required 
to resolve these features. Similarly, multi-frequency ra- 
dio observations can be used to probe scintillation ef- 
fects such as scatter-broadening. Higher-resolution spec- 
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troscopy will also be useful. For example, if the sin- 
gle emission line associated with the lens is [O II], giv- 
ing zi = 1.145, then moderate-resolution spectroscopy 
should split the doublet and allow an unambiguous line 
identification to be made. Finally, grid-based modeling 
of the resolved Einstein rin g flux may provide more ro - 
bust limits on substructure (jVegetti fc Koopm ans 2009), 
allowing us to better quantify the effects of mass pertur- 
bations in this system. 

5. CONCLUSIONS 

We use the Keck NIRC2 plus LGS AO to observe the 
gravitational lens system B0128+437, and we are able 
to confirm the presence of the lensing galaxy and clearly 
separate the lensing galaxy flux from the source galaxy 
flux for the first time. We take advantage of the high- 
quality AO imaging to present new science results for 
the B0128+437 lens system. These results can be sum- 
marized as follows: 

• For the most probable redshift zi = 1.145, we 
measure a total K-band luminosity log 10 (Lif) = 
(11.07 ± 0.05) h^Q L©, an aperture luminosity 
within the Einstein ring radius log 10 (L.K-,Ein) = 
(10.61 ± 0.05) L e , and a K-band M/L ratio 
M/L K ,Em = (0.68 ± 0.05) h 70 M o /L for the lens- 
ing galaxy. 

• We measure Kp flux ratios of Jb/Ia = 0.17 ±0.03, 
fc/fA = 0.42±0.01, and f D /f A = 0.30±0.03. The 
/b/ f A value in the near-IR is significantly smaller 
than both its expected value /b/Za ~ 1 and its 
measured radio values. This could be due to the 
presence of massive substructure, propagation ef- 
fects, or extragalactic microlensing. We find that 
none of these scenarios alone seems able to fully 
explain all aspects of the anomaly, implying that 
multiple phenomena are responsible. 

• We find that there is no evidence for luminous sub- 
structure above a point-source limiting magnitude 
Kp = 25.5, corresponding to My = —17 in the 
B0128+437 field. 

Looking at the wealth of new information obtained 
from B0128+437, it is clear that AO imaging is an 
important technique for studying gravitational lenses 
in new ways. With future surveys such as the Joint 
Dark Enery Mission (JDEM) and Large Synoptic Sur- 
yey Telescope (LSST ) expected to find ^10000 lenses 
([Marshall et al.l)2005[ ). of which an estimated 10% will be 
of the small-separat i on, low- mass lensing galaxy variet y 
(jTurner et all Il98i lOrban de Xivrv fc Marshall l200ll . 
AO imaging of lenses will also play an important role in 
analyzing a large subset of future data. 
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ing AO observing and NIRC2 data reduction and Leon 
Koopmans, John McKean, Phil Marshall, and Simona 
Vegetti for useful discussions. DJL, CDF, and MWA ac- 
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